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1.  INTRODUCTION 


The  objective  of  this  work  is  an  improved  model  for  high  frequency 
acoustic  backscatter  from  shallow  water  sediment  at  shallow  grazing  angles. 
Initial  development  of  a  model  was  accomplished  by  Boyle  and  Chotiros  in 
1992  (ARL-TR-92-12).^  This  work  included  studies  of  a  Biot  model  for  acoustic 
transmission  and  reflection,  originally  developed  by  Stern,  Bedford,  and 
Millwater.2  Also  included  were  the  results  of  a  comprehensive  compilation  of 
shallow  grazing  angle  backscatter  data.  The  compilation  suggested,  for  some 
sediments,  the  possible  importance  of  a  new  hypothetical  mechanism  for 
backscatter;  resonant  scattering  from  trapped  bubbles.  A  preliminary  model  for 
acoustic  backscatter  from  a  spatial  distribution  of  trapped  bubbles  in  sediment 
was  then  presented. 

This  report  contains  improvements  to  the  bubble  backscatter  model, 
including:  (1)  a  method  for  modeling  spherical  acoustic  sources  in  the  water 
column,  (2)  incorporation  of  the  Biot  model  in  describing  acoustic  propagation 
through  the  sediment,  and  (3)  an  improved  estimate  of  the  trapped  bubble  size 
distribution  function.  This  new  model  begins  with  a  calculation  of  the  pore 
pressure  in  the  sediment,  given  a  spherical  acoustic  source  In  the  water 
column.  This  procedure  is  discussed  in  Section  2.  Section  3  then  describes 
how  the  backscattering  strength  of  the  sediment  is  computed.  In  Section  4,  the 
resulting  theoretical  model  is  compared  with  experiment.  A  discussion  of  the 
results  follows  in  Section  5. 
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2.  PORE  PRESSURE 


In  this  report  we  consider  the  pressure  field  generated  by  a  point  source 
in  the  water  column  above  a  sediment  boundary,  illustrated  in  Fig.  2.1.  Of 
interest  is  the  fluid  pressure  in  the  interstitial  pores.  This  problem  is  treated  by 
decomposing  the  spherical  wavefield  into  a  spectrum  of  plane  waves.  The 
technique  follows  the  treatments  of  Brekhovskikh^  and  Westwood.^  A  Biot 
penetration  model  is  used  to  calculate  the  transmission  coefficient  for  each 
plane  wave.  The  pore  pressure  is  then  obtained  by  summing  the  contributions 
of  ail  plane  waves. 


2.1  PLANE  WAVE  DECOMPOSITION 


The  linear  acoustic  wave  equation  in  a  homogeneous  medium  with 
sound  speed  c  is 


2 

Vp  = 


2 

a  p 


2  ^2 

c  at 


(2.1) 


The  solution  to  the  above  equation  that  describes  the  pressure  field  caused  by 
a  spherical  source  is 


pKkR-toi) 

P  =  Po^-R—  »  (2-2) 

where  R  is  the  distance  from  the  source.  Po  the  pressure  amplitude  at  unit 
distance,  t  the  time,  cd  the  angular  frequency,  and  k=ci)/c  the  acoustic  wave 
number,  if  the  source  is  placed  at  the  origin  of  a  rectangular  coordinate  system, 
the  distance  R  is  given  by  R=(x2+y2+z2)^^2  |f  suppress  the 

time  dependence  e*'"’^  and  confine  ourselves  to  the  z=:0  plane,  the  pressure 
field  can  be  written  as 


ikr 

P(x.y)*^ 


(2.3) 
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where  r*s(x®+y2)^^. 

The  pressure  field  can  now  be  expanded  in  a  double  Fourier  integral  in 
the  variables  x  and  y. 

pikr 

—  =  03;iA(kx,ky)e'^'^*'^''^dkxdky  (r^tO)  ,  (2.4) 

where  A(kx,ky)  is  given  by 

(27c)^A(kx.ky)  =  dkxdky  (2.5) 

The  above  integral  is  more  tractable  in  polar  coordinates: 


kx  =  q  cos  y 


ky  s  q  sin  y 

q  =  (kx2+ky2)i/2 


X  srcos^ 

y  srsin^ 
dx  dy  =  r  dr  d(^ 


In  this  coordinate  system  A(kx.ky)  is  given  by 

(2n)2A(kx.ky)  =  dr d4> 


irtk-qcosW-e)] 
i(k -9008(^-6)] 


rspe 

d<^ 

r=0 


(2.6) 


if  there  is  some  absorption,  even  if  it  is  vanishingly  small,  k  will  have  a  positive 
imaginary  part.  In  this  case  the  quantity  in  brackets  is  zero  when  r  =  infinity. 
Therefore 


6 


WA(M,)=r 

j  0 

ii-aco»(#-0)) 

ft 


(2.6) 


The  integral  in  Eq.  (2.6)  is  tabulated.^  Upon  substituting  into  Eq.  (2.6),  the 
following  expression  for  A(kx,ky)  is  obtained: 


A(kx.ky)  = 


^k^-kx'ky 


(2.7) 


Substitution  of  Eq.  (2.7)  into  Eq.  (2.4)  yields 


2  ^ 

l\  *  IV^  •  l\y 


e'^*'^y)dkxdky  .  (2.8) 


The  above  expression  for  the  pressure  field  in  the  z^O  plane  can  be 
generalized  to  describe  the  field  in  all  space.  This  is  done  by  multiplying  the 

Integrand  by  e  for  z>0  or  by  e  for  z<0. 


p(x,y.z)  =  £[ 


2ic>/k^-k5-ky 


(z^O)  (2.9) 


p(x.y,z)  = 


V-k5-k2 


ei^*'**>'y'^**)dkxdky  ,  (z^O)  (2.10) 


where 


kj=Vk^-kJ-ky 


(2.11) 
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This  “continuation"  is  justified  because  the  integrand  still  satisfies  the  wave 
equation  and  because  the  new  expressions  reduce  to  the  original  expression 
Eq.  (2.8)  for  the  z^O  plane  when  z  approaches  zero. 

The  integrals  in  Eq.  (2.9)  and  Eq.  (2.10)  are  over  the  entire  kxky  plane.  In 
order  to  span  this  plane,  according  to  Eq.  (2.11),  imaginary  values  of  kz  must  be 
included.  These  so-called  “inhomogeneous  waves”  attuenuate  exponentially 
with  spatial  dimension  z  and  must  be  included  to  complete  the  construction  of  a 
spherical  wavefield  with  plane  waves. 

In  order  to  evaluate  the  integrals,  Eq.  (2.9)  and  Eq,  (2.10),  a 
transformation  to  spherical  coordinates  Is  made,  as  shown  in  Fig  2.2: 

kx  =  k  sine  cos(|>  . 

ky  =  k  sine  sin<j>  ,  (2.12) 

kz  =  k  cos6 


In  the  transformation  the  elemental  area  on  the  kxky  plane  becomes 

dkxdky  =  (k  sine  d<j>)(k  cose  de)  (2. 1 3) 

In  order  to  cover  the  entire  plane,  (j)  goes  from  0  to  2k.  By  allowing  e  to  run  from 
0  to  jc/2,  a  disk  on  the  plane  of  radius  k  is  covered.  In  order  to  span  the  rest  of 
the  plane,  q  is  further  allowed  to  vary  between  jt/2  and  jt/2-|oo.  It  would  also 
have  been  possible  to  span  the  plane  by  letting  q  vary  between  k/2  and  7i/2+ioo, 
but  this  would  have  made  the  Inhomogeneous  waves  increase  in  amplitude 
with  increasing  z  instead  of  attenuating. 

In  the  new  coordinates,  Eq.  (2.9)  becomes 


ei(k.x+isy4k,z)sineck|,de  .  (z^O)  .  (2.14) 

j^=o 
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in  order  to  determine  the  pressure  field  below  the  sediment  interface,  as 
in  Fig.  2.3,  each  plane  wave  component  must  be  multiplied  by  an  appropriate 
plane  wave  transmission  coefficient.  The  resulting  pressure  is  given  by 


_ 

p(x,y.z)  =  J^P  f  e'^'‘*’‘^'''^**’‘^“**^W(e)sineckj>d0,  (z^O)  ,  (2.15) 

where  W  is  the  transmission  coefficient,  which  depends  on  the  angle  of 
incidence  9.  is  the  height  of  the  source,  Z2  is  the  depth  of  the  receiver,  and  kzs 
is  the  wavenumber  in  the  lower  medium.  It  can  be  shown  that  the  integral  over 
<|>  is  a  Bessel  function.  This  is  done  by  first  expressing  the  variables  x  and  y  in 
polar  coordinates: 


x=rcos(|», ,  y=rsin(|)j 


The  integral  over  <|>  can  now  be  written  as 


-L 


=  f ^  gikrsine{coa>,coa|H«in*,8in^)^j. 
^  gikr  8in0(co8^,  co8(4h^,)]jj^ 


=  2nJ({krsin6) 


(2.16) 


(2.17) 


By  combining  Eq.  (2.14)  and  Eq.  (2.16),  a  simplified  expression  for  the  pressure 
Is  obtained: 


p(x,y,z)  =  1 2  Jo(kr  sine)  W(e)  sine  dS  .  (2.1 8) 

Je=o 


2.2  INCORPORATION  OF  BIOT  MODEL 

In  order  to  solve  Eq.  (2.18)  for  the  pressure  below  the  sediment 
boundary,  the  transmission  coefficient  W(0)  and  vertical  wavenumber  kzs  in  the 
sediment  must  be  specified.  For  this  purpose  the  Biot  penetration  model 
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described  in  ARL-TR-92-12B  was  used.  This  model  was  originatly  developed 
Stem,  Bedford,  and  Millwater^  and  later  modified  to  treat  gassy  sediments.  This 
was  done  by  varying  the  complex  compressibility  and  density  of  the  pore  fluid. 
The  gas  contribution  to  these  parameters  was  obtained  using  a  theory 
developed  by  Hawkins^  involving  reflection  coefficients  from  bubble 
distributions. 

The  Biot  model  differs  from  traditional  ways  of  treating  acoustic  media  in 
that  it  allows  for  two  phases:  a  solid  skeletal  phase  through  which  a  fluid  phase 
is  allowed  to  flow  freely.  It  is  based  on  the  following  coupled  pair  of  equations 
of  motion: 


+  (H-2|i)V(V-  u)  -  CV(V-  w)  =  pii  -  pfW  (2. 1 9) 

CV(V-  u)  -  MV(V-  w)  =  p|tt  -  (cpf/p)w  -  {Fti/k)W  .  (2.20) 

where  u  and  w  are  the  fluid  and  solid  displacements,  p  is  the  fluid  viscosity,  pf 
and  p  are  the  mass  densities  of  fluid  and  saturated  sediment,  respectively,  k  is 
the  permeability,  and  F  is  a  constant  that  accounts  for  the  frequency 
dependence  of  the  drag  force  between  fluid  and  solid.  The  constants  H,  C,  and 
M  are  determined  by  the  compressibilities  of  the  fluid  and  solid  constituents. 

One  important  feature  of  the  Biot  model  is  that  it  predicts  two 
compressional  waves  in  the  medium,  in  addition  to  a  shear  wave.  For  the  faster 
of  these  waves,  the  fluid  moves  generally  in  phase  with  the  solid  part  of  the 
medium.  For  the  slow  wave  the  two  components  move  out  of  phase.  Usually 
the  slow  waves  can  be  neglected  and  single  phase  models  can  be  used  to 
model  ocean  sediments.  Under  certain  conditions,  however,  slow  waves 
appear  to  be  dominant.  Such  is  the  case  when  a  wavefront  is  incident  upon  a 
sediment  half  space  at  a  grazing  angle  less  than  the  critical  angle  for  the  fast 
wave. 


The  above  described  Biot  penetration  model  calculates,  for  any  complex 
angle  of  incidence  0,  the  pore  pressure  transmission  coefficient  W(0)  and 
vertical  wavenumber  kzs.  The  transmission  coefficient  W(e)  is  obtained  from 
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simultaneous  solution  of  equations  (51-54)  of  Stem  et  at.  The  fast  and  slow 
vertical  wavenumbers  kzs  are  given  by  Eq.  (59)  and  Eq.  (60),  respectively. 

Equation  (2.18)  is  then  calculated  numerically,  given  W(6)  and  kzt- 
Figure  2.4  is  a  plot  of  the  computed  pore  pressure  in  the  sediment,  given  the 
input  parameters  of  Table  2.1.  The  pressure  contributions  are  broken  down 
according  to  the  types  of  waves  involved.  Wavespeeds  of  1690  m/s  for  the  fast 
wave  and  1170  m/s  for  the  slow  wave  are  predicted,  with  attenuations  of  12.2 
and  41 .2  dB/m,  respectively.  It  is  clear  that  the  slow  wave  has  a  pressure 
amplitude  that  is  not  negligible  in  comparison  to  that  of  the  fast  wave, 
particularly  when  the  horizontal  distance  to  the  acoustic  source  is  large.  Hence 
the  effects  of  the  slow  wave  must  be  included  in  predicting  total  pressures  in  the 
sediment  as  well  as  backscattered  pressures  from  the  sediment. 
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Scale  -  dB 


TABLE  2.1 

PENETRATION  MODEL  INPUT  PARAMETERS 


Fluid  Density 

1000 

Fluid  Bulk  Modulus 

2.0x109 

Pa 

Porosity 

0.4 

2650 

HE*SS9ii 

Pore  Size 

1.0x10-5 

m 

Viscosity 

1.0x10-3 

kg/m-s 

Permeability 

1.0x10-10 

m2 

Grain  Bulk  Modulus 

9.0x1 09 

Pa 

Frame  Shear  Modulus 

2.61x107 

Pa 

Shear  Logarithmic  Decrement 

0.15 

Frame  Bulk  Modulus 

4.36x1 09 

Pa 

Bulk  Logarithmic  Decrement 

0.15 

Virtual  Mass  Parameter 

1.25 

Gas  Bulk  Modulus 

2.48x105 

Pa 

Bubble/Pore  Radius  Ratio 

0.06 

Gas  Content 

0 
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3.  SEDIMENT  BACKSCATTERINQ  STRENGTH 

3.1  CALCULATION  OF  BACKSCATTERINQ  STRENGTH  BY 
RECIPROCITY 


Consider  the  situation  illustrated  in  Fig.  3.1.  In  Appendix  A  an 
approximation  for  the  backscattered  pressure  amplitude  P2  at  A  from  a  scatterer 
at  B  is  derived: 


P2  = 


(3.1) 


where  a  is  the  effective  cross  section  of  the  scatterer,  defined  as  the  area  which, 
when  multiplied  by  the  incident  intensity,  equals  the  total  scattered  power.  The 
incident  pressure  amplitude  generated  by  the  source  at  a  distance  of  1  m  in  the 
water  column  is  po.  and  pi  is  the  pressure  generated  at  the  scatterer,  located 
atB. 


In  this  analysis,  the  sediment  is  a  medium  with  many  scatterers 
throughout  its  volume.  A  mean  scattering  cross  section  per  unit  volume  Obv  can 
be  computed,  given  the  scattering  cross  sections  and  population  densities  of 
the  contained  scatterers.  If  a  projector  at  A  insonifies  a  volume  element  dxdydz 
at  B,  the  returned  pressure  at  A  is 


Obv dxdydz  p, 
4ic  Po 


(3.2) 


The  scattering  can  be  treated  as  a  function  of  location  along  the  water- 
sediment  interface,  as  in  Fig.  3.2.  in  this  case,  the  backscattered  intensity  from  a 
surface  element  dxdy  is  the  sum  of  all  backscatter  contributions  from  the 
sediment  volume  below  the  surface  element  dxdy.  The  total  mean  square 
backscattered  pressure  per  unit  surface  area  is  given  by 
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(3.3) 


Figure  3.3  illustrates  the  computation  of  the  backscattering  strength  of  a 
surface  element.  The  backscattering  strength  is 


BSs  lOlog 


(3.4) 


where  ps  is  the  root  mean  square  scattered  pressure  at  unit  radius  from  the 
surface  element  and  po‘  is  the  incident  pressure  at  the  surface  element.  The 
incident  pressure  po‘  is 


Po**Po/R  .  (3.5) 

where  R  is  the  distance  from  the  source  to  the  surface  element.  The  scattered 
pressure  at  unit  radius  is  obtained  from  the  scattered  pressure  Pbs  at  the  source, 

Pb  *  Pbs  R  (3.6) 

By  combining  Eqs.  (3.3)  -  (3.6),  the  backscattering  strength  Is  expressible  as  a 
function  of  obv.  Po,  and  pi. 


3.2  DETERMINATION  OF  SEDIMENT  VOLUME  SCATTERING 
CROSS  SECTION  Obv 

In  order  to  find  the  backscattering  strength  via  Eq.  (3.7),  it  is  necessary  to 
calculate  Obv.  the  backscattering  cross  section  per  unit  volume  of  sediment. 
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Since  the  hypothetical  mechanism  for  scattering  in  this  model  is  bubble 
resonant  backscatter,  obv  can  be  expressed  as  follows: 


Obv  =  J[*%(rb)c(rb)drb 


(3.8) 


where  fb(rb)  is  the  spatial  distribution  function  of  bubbles,  defined  as  the  number 
of  bubbles  per  unit  volume  per  unit  radius  increment  about  the  radius  rb.  and 
o(rb)  is  the  scattering  cross  section  of  a  single  bubble  of  radius  rb.  Expressions 
for  o(rb)  and  fb(rb)  are  developed  in  subsections  3.2.1  and  3.2.2. 

In  developing  the  model,  three  assumptions  are  made: 

(1 )  A  bubble's  size  is  regulated  by  the  pore  in  which  it  is  found;  i.e..  the 
ratio  of  bubble  radius  to  pore  radius  is  constant.  This  assumption  is 
necessary  to  permit  an  estimate  of  the  bubble  size  distribution  in  the 
absence  of  other  information. 

(2)  The  possibility  of  multiple  scattering  between  bubbles  is  neglected. 
Justifcation  for  this  assumption  is  included  in  Appendix  B. 

(3)  The  bubbles  act  as  though  they  were  spherical,  isolated,  and 
surrounded  by  water.  The  effects  of  the  interaction  between 
bubbles  and  between  bubbles  and  the  walls  of  the  pores  are 
neglected. 

There  is  some  experimental  and  theoretical  support  for  this 
approximation.  Howkinss  predicted  and  measured  the  resonance 
frequencies  of  bubbles  near  rigid  boundaries.  Two  configurations 
were  studied,  one  involving  a  bubble  near  a  semi-infinite  rigid  plane 
and  another  involving  a  bubble  contained  inside  a  rigid  cylindrical 
cup.  In  neither  case  was  there  a  dramatic  influence  on  the  bubble's 
resonance  frequency.  Resonance  frequencies  decreased  by  15% 
for  the  plane  boundary  and  50%  for  the  cylindrical  cup.  A  bubble 
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trapped  between  sand  grains  might  act  similarly,  being  more 
confined  than  a  free  bubble  near  a  semi-infinite  plane  but  less 
confined  than  one  inside  a  cylindrical  cup. 

3.2.1  Estimate  of  Single  Bubble  Beckecattering  Cross  Section 


Equation  (3.8)  requires  that  a(rb).  the  single  bubble  backscattering  cross 
section,  be  specified.  The  scattering  cross  section  for  a  single  spherical  bubble 
of  radius  rb  insonified  by  a  signal  of  frequency  f  is  given  by  Wildt;^° 


<Jr(rb)  = 


(3.9) 


5  is  a  damping  constant  and  fr  is  the  resonance  frequency,  given  by  Minneart^^ 


1 

27crb 


(3.10) 


where  y,  P.  and  p  are  the  ratio  of  specific  heats,  pressure,  and  density, 
respectively,  of  the  medium  inside  the  bubble.  Medwin^^  showed  that,  for 
resonant  bubbles,  rbfr  is  approximately  3  and  5  is  approximately  0.1.  It  is 
convenient  to  work  in  terms  of  the  variable  v  =  In(rb)  instead  of  rb.  In  this 
notation  the  single  bubble  scattering  cross  section  is 


Ov(v)  = 


(3.11) 


3.2.2  Estimate  of  Bubble  Size  Distribution 

Equation  (3.8)  requires  that  the  bubble  size  distribution  function  fb(rb)  be 
specified.  Our  approach  is  to  first  obtain  a  pore  volume  distribution  function.  A 
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bubble  distribution  function  is  then  found  by  assuming  that  bubbles  always 
occupy  a  constant  fraction  of  the  available  pore  space,  in  accordance  with 
assumption  (1)  stated  previously. 

Pore  Size  Distribution  Function 


Finney  and  Wallace^ ^  investigated  the  formation  of  pores  between 
randomly  packed  identical  spheres.  Figure  3.4  is  a  histogram  of  the  distribution 
of  interstitial  pore  radii  for  an  aggregate  of  spheres.  It  can  be  approximated  as  a 
bimodal  distribution  function  p(x,y)  where  x  and  y  are  the  natural  logarithms  of 
the  grain  and  pore  radii,  respectively.  For  any  value  of  x,  p(x.y)  is  approximated 
by  a  sum  of  two  normal  distributions. 


p(x,y)  =  A 


^bi 


(y-x-a,)*/2b^ 


+  B 


(y-x-ai)*/2»4 


.  (3.12) 


where  a,  is  the  mean  and  bj  is  the  standard  deviation  of  the  ith  peak. 

Sediments  have  a  measurable  particle  size  distribution  function  fx(x). 
where  x  =  log(grain  radius).  A  pore  size  distribution  function  fp(y)  is  estimated 
by  convolving  this  function  with  p(x,y): 

tp(y)  =  fx(x)*p(x.y)  .  (3.13) 

In  writing  Eq.  (3.13)  we  are  neglecting  the  possibility  that  particles  of 
varying  sizes  pack  differently  than  those  of  a  single  size.  For  example,  small 
particles  may  occupy  the  pores  created  between  larger  particles.  An 
assumption  is  also  mt^de  that  irregular  sediment  grains  pack  like  perfect 
spheres.  The  justification  behind  proceeding  with  these  assumptions  comes 
from  the  measured  porosities  of  sandy  sediments.  Porosities  between  0.36  and 
0.4  are  common  for  sandy  sediments.  This  compares  closely  with  a  porostiy  of 
0.363  for  randomly  packed  identical  spheres.^^ 
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Bubble  Size  Distribution  Function 


Assumption  (1)  states  that  a  bubble's  radius  is  governed  by  the  size  of 
the  surrounding  pore.  This  means  that  the  bubble  size  distribution  fb(v),  where  v 
is  the  natural  log  of  the  bubble  radius,  will  reflect  the  pore  size  distribution  fp(v); 

4(z)  =  ^(v-ln(Rbp))  .  (3.14) 

Rbp  is  a  constant  equal  to  the  ratio  of  a  bubble's  radius  to  that  of  its  surrounding 
pore.  It  is  assumed  that  each  pore  may  hold  only  one  gas  bubble,  and  that  the 
fraction  of  pores  that  indeed  have  gas  bubbles  is  In  the  absence  of  any 
further  information,  it  will  be  assumed  that  ^  is  a  constant  independent  of  pore 
size.  The  gas  fraction  ^  is  obtainable  by  dividing  the  total  volume  of  bubbles  by 
the  porosity: 


f  4(v)  dv 

C  =  - P -  •  (3.15) 

In  this  analysis  Rbp  and  C  are  treated  as  free  parameters  and  are  varied  to  fit  the 
model  to  experimental  data. 


4.  COMPARISON  OF  THEORY  WITH  EXPERIMENT 


in  Fig.  4.1  comparisons  are  made  between  the  predictions  of  the 
backscatter  model  and  experimental  measurements  of  the  backscattering 
strength  taken  over  three  selected  sites.  The  input  parameters  for  each  site  are 
listed  in  Table  4.1.  In  each  case  a  best  fit  was  obtained  by  varying  the  gas 
fraction  C  and  the  bubble>tO‘pore  radius  ratio  Rbp.  The  resulting  gas  fractions 
are  annotated  and  varied  between  1  x  10^  and  3  x  10*^. 

The  model  predicts  a  peak  in  the  backscattering  strength  at  a  frequency 
that  corresponds  to  resonance  for  the  bubbles  in  the  distribution.  The  width  of 
this  peak  will  depend  on  the  width  of  the  bubble  size  distribution.  The 
backscattering  strength  will  decrease  with  increasing  frequencies  that  are 
greater  than  resonance,  in  this  respect  the  model  differs  from  other  current 
models,  which  always  predict  an  increase  in  backscattering  strength  with 
increasing  frequency. 

At  present  the  model  is  still  being  developed  and  these  results  must  be 
regarded  as  preliminary.  An  impoved  estimate  of  the  compressibility  of  the 
gassy  fluid  phase  is  being  developed.  An  estimate  of  the  acoustic  impedance 
of  the  fluid  constrained  in  the  sediment  pores  is  also  being  incorporated.  These 
changes  to  the  model  will  probably  change  the  amplitudes  of  the  predicted 
backscattering  strengths  in  Fig.  4.1.  The  overall  features,  i.e.,  the  broad 
maximum  in  backscattering  strength  at  bubble  resonance  and  the  model's 
sensitivity  to  very  small  gas  fractions,  are  not  expected  to  change. 
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Kings  Bay.  1985 

Gas  Fractioa . 1.0x1  O'* 

Bubble/Pore  Radius  Ratio  (Rbp)...0.06 


0*3  10-2  10'^ 

Normalized  Grain  Size 


O'*  10'2  10'^ 

Normalized  Grain  Size 


Panama  City,  1986 

Gas  Fraction . 3x10'* 

Rbp . 0.1 


Jacksonville.  1 989 

Gas  Fraction . 2x10*^ 

Rbp . 0.01 


0-3  10'2  10'’ 

Normalized  Grain  Size 
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TABLE  4.1 

BACKSCATTER  MODEL  INPUT  PARAMETERS 


Site 

Kings  Bay 

1966 

Panama  City 
1986 

Jacksonville 

1989 

Ruid  Density 

(kg^m3) 

1000 

1000 

1020 

Ruid  Bulk 

(Pa) 

2.25x1 09 

2.25x10® 

2.406x10® 

Modulus 

Porosity 

0.3633 

0.3906 

0.3793 

Grain  Density 

{kg/m3) 

2656 

2650 

2650 

Grain  Radius 

(m) 

3.896x10*4 

1.593x10*4 

5.586x10*4 

Viscosity 

(kg/m-s) 

1.0x10*3 

1.0x10*3 

1.0x10*3 

Permeability 

(m2) 

1.0x10*10 

1.0x10*10 

1.0x10*10 

Grain  Bulk 

(Pa) 

9.0x1 09 

9.0x10® 

1.06x10l0 

Modulus 

Frame  Shear 

(Pa) 

2.61  x107 

2.61x107 

2.57x1 07 

Modulus 

Shear  Log 

0.15 

0.15 

0.15 

Decrement 
Frame  Bulk 

(Pa) 

4.36x10® 

4.36x10® 

4.36x10® 

Modulus 

Bulk  Log 

0.15 

0.15 

0.15 

Decrement 
Virtual  Mass 

1.25 

1.25 

1.25 

Parameter 

Gas  Bulk 

(Pa) 

2.48x1 05 

2.48x1 05 

2.48x1 05 

Modulus 

Bubble/Pore 

0.06 

0.1 

0.01 

Radius  Ratio 
Gas  Content 

1.0x10*5 

3x10*6 

2x10*6 
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5.  CONCLUSION 


A  heuristic  modei  for  high  frequertcy  acoustic  backscatter  from  sandy 
sediments  based  on  bubble  resonance  behavior  has  been  developed.  This 
modei  is  consistent  with  the  broad  backscattering  strength  maximum  observed 
in  some  of  the  experimental  data  and  demonstrates  that  a  trapped  bubble 
backscattering  mechanism  may  be  significant  in  high  frequency  acoustic 
backscatter  from  sediments. 

The  analysis  demonstrates  that  gas  bubbles  are  a  likely  cause  of 
observed  backscatter  from  sediments.  Very  small  gas  fractions  of  between 
1x10'^  and  3x10*6  are  necessary.  Since  there  is  no  direct  way  to  measure  gas 
fractions  this  small,  an  independent  verification  of  this  analysis  is  required, 
which  is  difficult. 

The  assumption  that  the  bubble  size  is  proportional  to  the  pore  size  was 
made  to  allow  a  first  estimate  of  the  bubble  size  distribution.  Further  research  is 
required  to  determine  a  more  precise  relationship  between  pore  geometry  and 
bubble  size. 

The  effects  of  Biot  slow  and  fast  compressional  waves  have  been 
included.  The  results  show  that  at  very  shallow  grazing  angles  the  slow  wave's 
amplitude  is  significant  in  comparison  to  that  of  the  fast  wave.  This  suggests 
that  a  complete  description  of  shallow  grazing  angle  penetration  must  include 
the  effects  of  the  slow  wave.  For  this  reason  the  Biot  theory  appears  to  be  most 
appropriate  in  acoustic  modeling  in  this  regime. 

At  present,  work  is  in  progress  in  improving  the  backscatter  model.  An 
expression  for  the  backscattering  cross  section  for  trapped  bubbles  is  being 
developed.  This  will  allow  elimination  of  the  assumption  that  trapped  bubbles 
scatter  as  though  they  existed  in  an  unbounded  medium.  Further  improvements 
will  include  studies  of  grain  packing  structures  and  their  relationships  to  trapped 
bubble  sizes  and  geometries. 
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APPENDIX  A 

RECIPROCITY  PRINCIPLE 
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Consider  the  situation  illustrated  in  Fig.  A.1.  According  to  the  reciprocity 
principle,  the  sound  pressure  at  B  when  a  source  is  placed  at  A  is  the  same  as 
that  produced  at  A  when  the  same  source  is  placed  at  B.  In  the  case  where 
sound  from  a  source  at  A  encounters  a  scatterer  at  B.  the  scatterer  can  be 
regarded  as  another  source  whose  sound  pressure  can  be  measured  back  at  A. 
In  this  case,  a  signal  propagating  along  the  return  path  experiences  the  same 
loss  that  the  original  incident  signal  did  on  its  way  to  B.  Let  a  cw  source  at  A 
cause  a  complex  pressure  pi  to  exist  at  B.  The  scatterer  at  B  will  return  a 
complex  pressure  P2  at  A,  given  by: 


where  piC  is  the  pressure  generated  by  the  scatterer,  measured  at  a  unit 
distance  from  the  scatterer.  The  incident  pressure  generated  by  the  source  at  a 
distance  of  1  m  in  the  water  column  is  po.  C  =  is  a  complex  factor  related  to 
the  scattering  cross  section.  Let  the  time  dependent  complex  pressures  po,  pi, 
p2  be  expressed  as 


Po  ®  Po® 


Pi  ®  Pi® 


j(6  i-ont) 


(A.2) 


P2  ®  P2® 


i(6ra)t) 


* 


where  pi,  P2.  and  p3  are  moduli  and  61.62.  and  av5  63  are  phases  of  the 
complex  pressures  po,  Pi,  and  P2.  Upon  substitution  into  Eq.  (A.1). 


P2® 


i(6fO>t) 


=  Pi® 


i(6  j-wt) 


re'1 


Pie 


i(6,-ci>t) 


iPo® 


i(0(ra>t) 


(A.3) 


Upon  taking  a  time  average  of  Eq.  (A.3), 


33 


Projector 

Backscattered  Pressure 
Amplitude  =  P2 


Scatterer 
Incident 
Pressure 
Amplitude  «  pi 


RGURE  A.1 

BACKSCATTER  AT  A  FROM  SINGLE  SCATTERER  AT  B 


ARL-UT 

AS-93-473 

NPC/RFO 

8-93AkC 


34 


The  factor  r  is  related  to  the  scattering  cross  section  o  as  follows: 


r  = 


(A.5) 


where  o  is  defined  as  the  area  which,  when  multiplied  by  the  incident  intensity, 
equals  the  total  scattered  power.  Upon  substituting  Eq.  (A.5)  into  Eq.  (A.4)  and 
noting  that  the  time  averages  of  the  bracketed  quantities  in  Eq.  (A.4)  are  unity. 


P2 


Y  4n  Po 


(A.6) 


The  situation  in  Fig.  3.1  differs  from  that  in  Fig.  A.1  in  that  the  scatterer  exists 
under  a  fluid-sediment  boundary.  The  above  derivation  leading  up  to  Eq.  (A.6) 
rests  on  the  assumption  that  both  scatterer  and  projector  exist  in  the  same 
medium.  For  the  case  in  which  the  scatterer  is  below  a  sediment  boundary,  an 
impedance  mismatch  between  fluid  and  sediment  must  be  accounted  for.  In 
this  analysis,  such  a  mismatch  is  neglected.  As  a  result,  Eq.  (A.6)  can  only  be 
considered  an  approximation.  The  approximation  should  be  reasonable  since 
the  trapped  bubbles  that  act  as  scatterers  exist  in  the  trapped  fluid  within  the 
sediment  pores.  Effectively,  the  assumption  is  being  made  that  the  interstitial 
fluid  has  the  same  acoustic  impedance  as  the  fluid  above  the  sediment. 


This  pags  intsntionally  left  blank. 
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APPENDIX  B 

JUSTIFICATION  OF  SINGLE  SCATTERING  MODEL 
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Attenuation  in  sediment  is  caused  by  absorption  and  scattering  of  an 
incident  acoustic  wave.  The  single  scattering  approximation  is  valid  when  the 
absorption  contribution  is  large  in  comparison  with  the  scattering  contribution. 
In  practical  terms  this  criterion  is  considered  satisfied  when  the  single  scattering 
albedo^  is  less  than  0.2.  The  albedo  a  is  defined  as 


a  = 


o 


(B.1) 


where  Os  is  the  scattering  cross  section  per  unit  volume  and  Oe  is  the  extinction 
cross  section  per  unit  volume.  Each  bubble  in  the  sediment  will  act  as  a 
scatterer  with  scattering  cross  section  Os  given  by  Eq.  (3.9); 


Os(rb)  = 


(B.2) 


The  total  bubble  scattering  cross  section  per  unit  volume  Os  is  this  expression 
integrated  over  the  distribution  of  bubble  radii: 


i=f  a^rb)f(rb)drb 

•/  0 


(B.3) 


The  extinction  cross  section  Oe  associated  with  the  same  bubble  is^’' 


<yc(rb)  = 


where  Kr  is  the  acoustic  wavenumber  at  resonance, 


(B.4) 


(B.5) 
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The  total  extinction  cross  section  per  unit  volume  is  the  integral  of  Eq.  (B.4)  over 
the  bubble  size  distribution. 


^=)  <y,(*‘b)f(rb)drb  , 

Jo 


(B.6) 


Upon  inserting  expressions  (B.3)  and  (B.6)  into  expression  (B.1)  above, 
the  single  scattering  albedo  is  given  by 


f  a,(rb)  f(rb)  d  Fb 
Jo 

f  o/Tij)  f(rb)  d  Fb 
Jo 


(B.7) 


The  integral  in  the  numerator  can  be  expressed  as  a  sum  of  resonant  and 
nonresonant  contributions.  The  nonresonant  part  is  negligible  except  in  the 
case  of  very  large  bubbles.  The  remaining  resonant  contribution  includes  a 
delta  function  about  the  resonance  radius.  The  value  of  this  integral  can 
therefore  be  approximated  by  simply  evaluating  the  integrand  at  the  resonance 
radius. 


Computation  of  the  albedo  via  expression  (B.7)  is  nontrivial.  In 
determining  whether  a  single  scatter  model  applies,  however,  it  suffices  simply 
to  determine  if  the  albedo  is  less  than  0.2.  The  numerator  can  be  estimated 
accurately  by  neglecting  the  non  resonant  contribution  as  described  above.  If 
the  denominator  is  approximated  similarly,  i.e.,  by  evaluating  its  integrand  at  the 
resonance  radius,  any  error  incurred  will  make  the  denominator  smaller.  An 
approximate  albedo,  which  can  be  defined  as 


~  <y.(ar) 

a  = -  , 

cyjtar) 


(B.8) 


will  therefore  be  greater  than  or  equal  to  the  actual  albedo  given  in  expression 
(B.7).  Therefore,  if  the  approximate  albedo  a  is  less  than  0.2,  a  single  scattering 
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model  is  sufficient.  Upon  combining  Eq.  (B.2).  Eq.  (B.4),  and  Eq.  (B.8). 
becomes 


- 


(B.9) 


Medwin^^  showed  that,  for  resonant  bubbles,  frdr  is  approximately  3  and 
S  is  approximately  0.1.  If  these  values  are  used,  a  is  evaluated  at  0.126.  A 
single  scattering  model  is  therefore  appropriate. 
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